On the third order structure function for rotating 3D homogeneous turbulent flow 
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A form for the two-point third order structure function has been calculated for three dimensional 
homogeneous incompressible slowly rotating turbulent fluid. It has been argued that it may pos- 
sibly hint at the initiation of the phenomenon of two-dimensionalisation of the 3D incompressible 
turbulence owing to rotation. 

PACS numbers: 47.27.i, 47.27.Jv 



I. INTRODUCTION 

With possible realistic applications to flows in oceans 
and atmospheres, rotating turbulence attracts interests 
of oceanographers, geophysicists, mathematicians and 
physicists alike. Rotation seems to be serving as a bridge 
between 2D, quasi-2D and 3D turbulences, a fact well- 
established with the discovery of two-dimensionalisation 
of 3D turbulence due to rotation. 

In the steady non-turbulent flow, for low Rossby num- 
ber (Ro = U /2L£l) and high Reynolds number (Re = 
ULjv), Taylor-Proudman theorem[l] argues that rota- 
tion two-dimensionalises the flow. This argument is of- 
ten mistakenly extended to turbulent flows to explain the 
rotation induced two-dimensionalisation arising therein. 
The two-dimensionalisation of the 3D turbulent flow in 
presence of rotation has begun to be understood as a 
subtle non-linear effect which is distinctly different from 
Taylor-Proudman effect. 

Cambon et al. [2] showed that in the presence of rota- 
tion, the transfer of energy from small to high wavenum- 
bers is inhibited; at the same time, the strong angular 
dependence of this effect leads to a draining of the spec- 
tral energy from the parallel to the normal wave vectors 
(w.r.t. the rotation axis) showing a trend towards two- 
dimensionalisation. 

Waleffe[3] used helical decomposition of the velocity 
field to study the nature of triad interactions in homo- 
geneous turbulence and coupling it with the instability 
assumption predicted a transfer of energy toward wave 
vectors perpendicular to the rotation axis under rapid 
rotation. The helical decomposition turns out to be very 
handy to deal with rapidly rotating turbulent flow. In 
that case the linear eigensolutions of the problem, the so- 
called inertial waves, have the structure of helical modes. 
The assumption about the triadic transfers, coupled with 
resonance condition for non-linear interaction between 
inertial waves, show that there will be a tendency to- 



ward non-linear two-dimensionalisation of the flow. 

Simulations by Smith et a/. [4] speak volumes for the 
two-dimensionalisation effect. They showed the coexis- 
tence of inverse cascade (a typical feature of 2D turbu- 
lence) and forward cascade in forced rotating turbulence 
within a periodic box of small aspect ratio. In the simu- 
lations, the ratio of the mean rates of energy dissipated 
to the energy injected decreased almost linearly, for Ro 
less than a critical value, with decrease in Ro (increase 
in angular velocity |f2|). By the way, a very recent nu- 
merical study [5] shows similar transition from stratified 
to quasi-geostrophic turbulence, manifested by the emer- 
gence of an inverse cascade - a conclusion that agrees 
with that of Lindborg[6]. 

Although recent experiments by Baroud et al.[7, 8] and 
Morize et al.[9, 10] have shed some light on the two- 
dimensionalisation effect, the scaling of two-point statis- 
tics and energy spectrum in rotating turbulence remains 
a controversial topic. Zhou[ll] in analogy with MHD 
turbulence proposed an energy spectrum E(k) ~ fc~ 2 for 
rapidly rotating 3D turbulent fluid and this does seem 
to be validated by some experiments [7, 8] and numerical 
simulations [1 2-1 5]. But some experiments [9] do not tally 
with this proposed spectrum. They predict steeper than 
k~ 2 spectrum and this again seem to be drawing some 
support from numerical rcsults[I6, 17] and analytical re- 
sults found using wave turbulence theory [18, 19]. 

To be precise, if one wishes angular velocity to become 
a relevant parameter in constructing the energy spectrum 
E(k), simple dimensional analysis would lead one to: 



E(k) oc rr 



(1) 
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where m is a real number, m should be restricted within 
the range 5/3 to 3 to keep the exponents of O and e in 
relation (1) positive. The two limits m = 5/3 and m = 3 
corresponds to isotropic homogeneous 3D turbulence and 
2D turbulence respectively. The spectrum due to Zhou - 
E(k) ~ k~ 2 - is due an intermediate value of m — 2. So, 
as far as the present state of the literature on rotating 
turbulence goes, two-dimensionalisation of 3D turbulence 
would mean the dominance of a spectrum which goes to- 
wards E(k) ~ k~ 3 and which may choose to settle at 
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E(k) ~ k~ 2 , an issue yet to be fully resolved. 

These all studies are for low Ro high Re limit while 
high Ro high Re limit has been rather less ventured in 
relation to two-dimensionalisation of turbulence. Now, 
the two-point third order correlation function (S3) in ho- 
mogeneous isotropic turbulence has a rare non-trivial ex- 
act result. In this paper we have studied the two-point 
third order correlation function (5*3) in this latter regime 
of high Ro high Re and went on to argue that a spectrum 
which goes as k~ 2 - 87 (i.e., in between k~ 2 and fc~ 3 ) exists 
although dominated by the 3D spectrum E(k) <~ k~ 5 / 3 
the reason for which of course is that we arc dealing with 
low rotation limit. This obviously hints at the initiation 
of the two-dimensionalisation effect with slow rotation. 

Also, if one goes by the procedure given in the book by 
Frisch[20] to derive the form of the correlation function 
in d-D turbulence with the assumption of forward energy 
cascade, one would land up on [21]: 



-1 3v 
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d(d + 2) 



el (2) 



where e is the mean rate of dissipation of energy per 
unit mass. This result is not quite true for the two- 
dimensional case since it gives for d = 2, S3 = — (3/2)el 
and not S3 = (3/2)el because the calculation doesn't 
take into account the conservation of enstrophy which 
causes the reverse cascade of energy [2 2]. It might be 
noted that 5*3 = (3/2)el for d — 2 is for the regime of 
scales larger than the forcing scale[23, 24]. If, using cal- 
culations of structure functions, in the limit of high Ro 
and high Re, one wishes to see whether a trend towards 
two-dimensionalisation of 3D homogeneous isotropic tur- 
bulence occurs or not, then basically one would have to 
check (a) if S3 = — (4/5)el at small scales for 3D turbu- 
lence shows a tilt towards S3 = (3/2)el at large scales 
for the 2D turbulence and (b) if the forward energy cas- 
cade is depleted at the smaller scales. As we shall show 
here, in the lowest order calculation this is what one may 
get, again hinting at the initiation of the effect of two- 
dimensionalisation of 3D turbulence owing to the small 
anisotropy induced by slow rotation. 



II. RELEVANT SCALES IN ROTATING 
TURBULENCE 



unit mass), Q (angular velocity) and Iq (integral scale 
which typically is the system-size). The three important 
time-scales involved in the system are: U <~ e -1 / 3 / 2 / 3 
(eddy-turnover time or circulation time for the eddy of 
scale I; I < lo), tn <~ and td ~ I 2 ' jv (diffusion time 
scale). It is well-known that a length scale In = \J (e/^ 3 ) 
is what responsible for the estimation of the anisotropy 
introduced by the rotation. The competition between 
the time-scales t\ and i<j gives rise to what is known as 
dissipation length scale Id, defined as Id = (v/e) 1 ^ and 
a similar competition between the time-sc ales td and tn 
allows us to define a length scale Ind = \/ (v fVi). Now, 
lets look at the typical scenario when Ro is moderate. 
The four vital length scales are typically arranged ac- 
cording to the order : l n > In > Ind > Id- Thus, the 
regime lo > I > In is the regime where effect of rotation 
is important and anisotropy reigns. The scales I £ (Iq, Id) 
may be considered to have isotropy, though to be precise, 
probably Id here should be replaced by Ind since rotation 
seems to be bringing the effect of viscosity to rather larger 
length scales. So, now what happens when the Ro is de- 
creased by increasing the angular velocity is interesting. 
Both the scales In and Ind rush towards the dissipation 
length scale, thereby increasing the anis otropic regime 
and at the angular velocity Q — fl a = \J (sjv) one has 
Iq = hid = Id and the turbulence is fully anisotropic. 

Strictly speaking, even a small rotation introduces 
anisotropy (however small) at all scales and the isotropic 
regime does have a degree of anisotropy in it as we shall 
see in this paper. In the fully anisotropic limit, i.e. for 
= O a , one expects full decoupling of the plane per- 
pendicular to the rotation axis from the direction of the 
rotation axis. However, even in the partially anisotropic 
limit (e.g. when we have slow rotation imparted on the 
turbulent fluid), l z should still be given a special status 
for being in the direction of the rotation axis, by which 
we mean that the structure functions should no longer 
depend on I but rather on l z and l± (where I 2 = I 2 + 1\ 
and |0| = Q z ). 

We shall see how this decoupling sets in, in the limit 
of low angular velocity and try to study in that very 
limit, the two-point third order structure function in 
the first approximation and see how the effect of two- 
dimensionalisation is all set to sneak in with the switch- 
ing on of rotation. 



Let us look in to the various length scales that have 
to be taken into consideration while talking about a ho- 
mogeneous rotating turbulence which basically satisfies 
following version of Navier-Stoke's equation: 

^ + (V.v) v = -Jvf - d x (n x xj 

-2fZ x v + vV 2 v + f (3) 

Various parameters to be considered are: v (kinematic 
viscosity) , e (finite mean rate of dissipation of energy per 



III. S 3 FOR SMALL ft 

Let us start with low fi-limit. With this statement we 
mean, as discussed in the previous section, ft <g; il a . So, 
the entire fluid may still be treated as isotropic but as 
rotation should play a role, we assume that < ViVjv' k > 
(where angular brackets mean ensemble average and Vi = 
Vi(x,t) is the z-th component of velocity and similarly, 
v[ = Vi(x+l, t)) should depend on f2 as well. f2 would take 
care of the mild anisotropy. Since, physically speaking, 
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S3 should not depend on which way the rotation axis is 
and since we are interested in low values of Q, we shall 
let < ViVjv' k > depend only on the terms quadratic in 
n and not bother about higher order terms in f2. As a 
result, we write the following most general tensorial form 
for < ViVjv' k >: 



bi 3 .k = < ViVjv k > 

= C(l)6 t3 l° k + D(l)(5 ik l° + 5 0k l°) + F(l)l?l]l° k 
+G(l)[(e imk l° + e jmk l?)l° m } + HiW^ll 

+K{i)(n i n k i° + n j n k i°) (4) 



where 1° is the i-th component of the unit vector along 
I. We have assumed that the coefficients are dependent 
only on I and it is the f2 which is taking care of the 
mild anisotropy which the turbulent fluid might have. 
We must accept that the assumption of letting coeffi- 
cients depend only on I is rather crude in the light of 
the complex forms that the two-point tensors in a fully 
anisotropic turbulence flow take[25]. The justification, 
and hence solace, for the assumption, however, can be 
drawn from the fact that very simple revealing results 
matching with recent experiments are arrived at in the 
long run. As we are considering incompressible fluid, we 
must have: 



dlb,. 







(5) 



which when applied to relation (4), yields relationships 
between various coefficients. Einstein summation con- 
vention has been extensively followed in these calcula- 
tions unless otherwise specified. Using relations (4) and 
(5), one lands up in the end on the following: 

Bijk = ((v'i - Vi)(v'j - Vj )(v k - v k )} 
= 2{hj,k + b jk a + b ki j) 
= -2(lC + C){5 l3 l° k + b ik l° + 5 3k l°) 
+6(IC> - C)l°l°l° k 

+ui(£i i n j i k + n t n k i° + n 3 n k i°) (6) 

Here, prime (""') denotes derivative w.r.t. I and J is a 
constant which, curiously enough, is of the same dimen- 
sion [L 2 T _1 ] as that of the kinematic viscosity. Now we 
can see that using the relation (6), two-point third order 
structure function (5*3) can be extracted from Bij k in the 
following way: 



S 3 (l) EE ((SV^I)) 3 ) 



^v(x + I) — . 



Ss(l) = ((H ~ ViKWj ~ Vj)l?W k - v k )l° k ]) 
Ss(l) = B ijk l°l°l k 

s 3 (i) = -12C + — -(n.i) 2 



(7) 



where we have used relation (6). One may define physical 
space energy flux (e(l)) as: 



e(l) = -^Vi. < \5v(l)\ 2 6v(f) > 



(8) 



=* e(l) = lC" + 7C' + ^+3,m 2 + ^r(n.l) 2 (9) 

To get relation (9), we have again made use of the rela- 
tion (6). The energy flux through the wave number K 
(JIk) for the isotropic homogeneous turbulence may be 
calculated to be: 



' dii^S.(i + ia,)f(!) (10) 



Now if one makes the standard assumption (often made 
made during the derivation of S3) that as Re — > 00, the 
mean energy dissipation per unit mass e(v) tends to a 
positive finite value (i.e., \im. v ^Qe(v) = e > 0), then 
linv^o rix = £ in the inertial regime. Therefore, in the 
inertial range, putting x = Kl, one has 



2 f c 
n Jo 



dx 



sin(a;) 



f (i) 



where, 



/(J) =/(0 = (l + ldi)e(l) 



ill) 



(12) 



For small I (large K), the integral in relation (11) yields 



f(l) « e 



(13) 



Now using relations (9), (12) and (13), we form a differ- 
ential equation which when solved keeping in mind that 
53 should not blow up at I — one gets following form 
for S3 in slowly rotating homogeneous turbulent fluid. 



Ss(l) 



1 9 

^Ji[n 2 



m k i° k ) 2 } 



(14) 



One may note from the relation (14) that how magically 
has brought up the anisotropic effects even for small 
f2 though for the entire calculation we followed the pro- 
cedure for the homogeneous isotropic turbulence. Thus, 
the form for S 3 is pausible. 

One may ask: does the effect of two-dimensionalisation 
shows up in the relation (14)? As one may note from the 
relation (14) this is quite a possibility but the only catch 
being that J should be positive, an issue which we have 
not been able to resolve. If J is positive, it means if we 
increase f2 the value of S3 would distort away from the 
usual — (4/5)d for the non-rotating case to more positive 
values. This apparently shows that the effective value of e 
is decreased depicting that the forward energy transfer is 
depleted which is in keeping with the discussion given in 
the last paragraph of the section (I) and hence the ten- 
dency of the rotating 3D turbulence to show the effect 
the two-dimensionalisation is being highlighted. That 
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the sign of J should be positive is a question remains to 
be addressed. 

By the way, the relation (14) also suggests that the 
coefficients in the tensorial form for bij_ k should have de- 
pendence on l z and l± separately effecting a mild decou- 
pling of directions. So taking hint from it, we proceed 
to rewrite bij ik for slowly rotating 3D turbulent fluid but 
now introducing anisotropy directly into the coefficients 
and not letting il take care of anisotropy explicitly. Of 
course, the coefficients will now depend on ft. 

For completely isotropic homogeneous turbulence, one 
would write following general form (relation (15)) for bij tk 
which is made up of Kronecker delta and components of 
the unit vectors l/\l\. 

b ijtk = C{l)5 ij l° k + D{l){5 ik l° + 5 jk ll) 

+F(l)l°l°l° k (15) 

The expression is symmetric in i and j and the coeffi- 
cients are dependent on I only. As discussed earlier, with 
rotation coming into effect, anisotropy comes into effect. 
If this effects in the possible decoupling (even if partial) 
of the direction along the rotation axis (which we shall 
take along the z-axis), then mathematically we may in- 
troduce this effect by modifying the form (15) of bij ik to 
the following: 

bij, k = C(l, l z , n)S l3 l° k + D(l, l z , fl)(S ik l] + 6 jk l°) 

+F(l,l z ,Cl)l°l°l° k (16) 

If one uses the incompressibility condition (relation (5)), 
one gets: 



D= l -(-C'-^)-C 

b = o 



(17) 

and D = (18) 

where dot represents the derivative w.r.t. l z and prime, 
as before, the derivative w.r.t. I. Using equation (17) in 
the equation (18), one land up on: 

Cl z + IC' + 3C* = 
=> C = Y^A n r n - 2 l n z 

n 

=> C ^ for n £ (-oo, -2] n [0, oo) 

=> C = D = F = Q (19) 

In arriving at the result (19), we have taken care of the 
fact that C can not be allowed to blow up for either for 
l z = or for I = 0. Thus the relation (16) vanishes 
trivially. So, we are left with the following choice: 

bij, k = C(l±,l z , Sl)Sijl% + D(l ± ,l z , tt){5 lk l° + 5 lk l°) 
+F(l ± ,l z ,Sl)l°l°l° k (20) 

Using relation (5) and relation (20), we arrive at following 
relationship between the coefficients: 

D = -—C--C-C (21) 
2 2 v ; 

F =P + lt 6+ (M-T) 6 -^- c (22 » 



Here tilde and dot define derivatives w.r.t. l± and l z 
respectively Proceeding monotonously as before we get 



Bijk — 2(bij,k + bjk,i + bhi,j) 

= -2(l ± C + l z C + d)(Sijll + 6 ik l° + 5 jk l°) 
+6Fl?l°l° k 

And hence, 



(23) 



S 3 = B ijk l?l?l° k = 6[F - {l ± C + l z C + C)} (24) 

The definition for the physical space energy flux (e(l)) 
has to be obviously modified to: 

< \5v(l)\ 2 8v(l) >= B iia l°J± + B iiz l° l f (25) 

t_L lz 

where a takes two values - x and y only. Now, using 
relations (8), (22), (23) and (25) and performing tedious 
algebra one gets: 



-1 



(mi + erji + siir z )c 



4(ii + p z y 

+ (3l 5 J z + 6l 3 J 3 z + 3l 4 J 2 z + 6l 2 J 4 z + 31 Jl 

+3l 6 z ) C + (3l 3 J 3 z + 6l±l 5 z + C 
+ (5l 5 ± + 6l 4 J z + 2il\l 2 z + 12l 2 J 3 z + 181 J 4 Z 

+6l 5 z ) C + {-7l\l z + U\l 2 z - l 2 J 3 z 
+231 Jl + 6l 5 z + mi 1 ^) t + (-12l 4 z 
-8l 3 J z - 20l 2 J 2 z + 361 J 3 Z + 181 A Z 
+8l~ 1 l 5 z ) C + (-13l 3 J z - 43l 2 J 2 z - 39l ± l 3 z 

-17/*) C + (-4/i - 8l\l z - \2lJl) C 

(26) 

The energy flux (11^-) through the wave number K for 
the homogeneous (not necessarily isotropic) turbulence 
may be shown to be: 



1 

2^2 



(27) 



Using cylindrical polar coordinates we reduce the relation 
(27) to: 
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l±dl±dl z <{ — - — -x 



l 2 dl ± l 2 dl z i 2 



e(l) (28) 



Now, we introduce the variables y = Kl± and z = Kl z 
in relation (28) to get: 



/"OO POO 

n K = - / 



dydz 



sin(y 2 + z 2 ) 2 



f(- ~) 



(29) 
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Now, let's probe small I behaviour. Because 

ir=-ooIyZ dydz[ S m(y 2 + z 2 ) 1/2 ]/{y 2 + * 2 ) = - 2 A we 

have 



7T 



(30) 



Obviously, s has the meaning of finite positive mean rate 
of dissipation of energy per unit mass. Using the expres- 
sions (26) and (30), we look for the l z = limit. One 
then has the result: 



d 
dl ± 



1 



(3*3.(5 + U^C - 12(7 - 4^) = - — 
8e, 



3l 4 ± C + Ul 6 ± C - 2l 2 ± C - 161 ± C = l ± 

IX 

/ 7-V97 T+VWf 

C= [A 1 + A 2 ll 1 + A 3 l ± e +A 4 l ± B 



+ 



el±_ 
~2t7 



(31) 



Relations (22), (24) and (31) together yield following ex- 
pression for 53: 



S: 



3\l z =0 



6 / 

-— £tj_ 

7T 



-eh 









K- 




7+V97 


12\l ± « 


7+V97 



1 + V97 



Al ± 



(32) 



where, A is a constant which for obvious reason depends 
on Q and e. Using dimensional arguments and introduc- 
ing a non-dimensional constant c, we may set 



. „ 1 + Vof lWg 

A = Cll 4 £ 12 



(33) 



From relations (32) and (33), we may write finally 

<S*3 j i_=0 



6 I + n/97 11-V97 i±vH 

-el± + cn—i—e—ti—L 6 (34) 

TT 



This (relation (34)) is the final form for two-point third 
order structure function in the plane whose normal is par- 
allel to the rotation axis for slowly rotating homogeneous 
3D turbulence. 



IV. ENERGY SPECTRUM FOR SMALL Q 

If we for the time being forget about the issue of 
anomalous scaling, then a mere inspection of the rela- 
tion (34) from the point of view of dimensional analysis 



would tell that in the directions perpendicular to the axis 
of rotation, there are two possible energy spectrums viz. 



and, 



E(k) 
E(k) 



fc-S 

16+V57 



(35) 
(36) 



which are respectively due to the first term and the sec- 
ond term in the R.H.S. of the relation (34). It is very 
interesting to note that the exponent of k in the relation 
(36), i.e. -(16 + V97)/9, equals -2.87 which is in be- 
tween —3 (for 2D turbulence) and —2 (for rapidly rotat- 
ing 3D turbulence as proposed by Zhou). Obviously, the 
spectrum (35) will be dominant compared to the spec- 
trum (36). But as the is increased (of course, remaining 
within a range so that the anisotropy is not strong enough 
to breakdown the arguments used to calculate the S3 of 
the relation (34)), the spectrum (36) becomes more and 
more prominent; thereby two-dimensionalisation of the 
3D homogeneous turbulent fluid is initiated which then 
carries over to high rotation regime as is being exten- 
sively studied. This signature of two-dimensionalisation 
is, of course, in agreement with what present literature 
on turbulence hails as the two-dimensionalisation of tur- 
bulence. 



V. CONCLUSION 

To conclude, we emphasis on the fact that the 
form of two point third order structure function in 
a slowly rotating homogeneous 3D turbulence can 
strongly hint towards the initiation of the effect of two- 
dimensionalisation of 3D turbulence. It barely needs to 
be mentioned that the relations are quite interesting and 
pertinent (at least within the approximations made in the 
calculations) - something which is worth getting in the 
literature of turbulence since exact relations are very few 
therein. So any theory developed in the limit of Ro — > 
and Re — ► 00, must satisfy the relation derived in this 
paper in the limit of low f2 or explicitly violate the as- 
sumptions made to arrive at the result; in this sense the 
relation (34) may prove to be of high importance. 
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